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Clase 1 



Generalidades 

Los campos electrico (E) y de induction magnetica (B) se introdujeron orig- 
inalmente a traves de la fuerza ejercida por cargas (q r ) o corrientes (J) sobre 
una carga de prueba (g): 

F q = qE <-> E 
dFj = Idtx B <-> B 

— * 

De acuerdo con esto, B se interpreta como fuerza por unidad de carga y B 
como fuerza por unidad de corriente. Sin embargo ambos campos tienen 
significado propio, independiente del tipo de fuente que los genera. 

Ahora bien, E y B no son los linicos campos importantes en la elec- 
trodinamica. En la mayoria de las sustancias: 

D = e E + P, 

H = —B-M . 
Mo 

— * 

D se conoce como desplazamiento electrico, H campo magnetico, P y M son, 
respectivamente, las polarizaciones electrica y magnetica (i. e., representan 
el pro medio macroscopico de dipolos electricos/magneticos en el material en 
presencia de campos); e = 8.85 x lO~ 12 C/Nm 2 , fi = Anx 1CT 7 H/m y e /io = 
c~ 2 . La conexion entre los vectores (P, E) y (M, H) esta determinada por las 
propiedades de cada sustancia. Para medios anisotropicos la aproximacion 
lineal en los campos es: 

Pi = e aikE k , 
Mi = K ik H k , 

con i, k — 1, 2, 3; a es el tensor de polarizabilidad y k el tensor de magneti- 
zacion. Entonces 

Di = €ikE k , 
Bi = fiikHk, 



\dF T 



Idl sin a 
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donde e ik = e (5 ik + a ik ), (x ik = Ho{8 ik + K ik ). Para medios isotropicos: 

P = e aE 

M = kH, e = e (l + a) 
D = eE, fj, = /i (l + k) 
B = HH. 

Una vez definidos estos vectores, podemos presentar las ecuaciones de 
Maxwell (1873), que son en electrodinamica lo que las leyes de Newton en 
mecanica. Las ecuaciones de Maxwell en forma diferencial son 

- dD 

Vxtf= J + -, (!) 

V • B = 0, (2) 

Vx*=-f. (3) 

V • D = p (4) 

(p es la densidad de carga y j la densidad de corriente). La forma integral 
de estas ecuaciones es 



B ■ MA = 

s 



( pdV. 
Jv 



D ■ MA 

De estas ultimas se obtienen las condiciones de frontera entre dos medios: 

(D 2 - D 1 ) ■ ni, 2 = a, (5) 

ni, 2 x (Ex - E 2 ) = 0, (6) 

(B 2 - B x ) ■ n 1)2 = 0, (7) 

ni, 2 x(H 2 -H 1 )= I (8) 



3 



donde 

N " Is ■ 

Electrostatica 

Estudiamos un problema de electrostatica si se satisfacen las condiciones 

• No hay dependencia temporal en los campos. 

• No existen cargas en movimiento. 

Con esto, las ecuaciones de Maxwell ([l] - §) se reducen a 

Vx£ = 0, (9) 
V-D = p. (10) 

En vista de (||), del calculo vectorial sabemos que 

E = -V$. 

De esta forma se introduce el potencial electrostatico ($). Considerando 
medios isotropicos (i. e., D = eE) la ecuacion (pi]) se reduce a 

V 2 $ = - P - 

e 

que se conoce como ecuacion de Poisson (en ausencia de cargas se obtiene la 
ecuacion de Laplace). 

Por otra parte, las condiciones de frontera (|5| - H) se reducen a 

$1 = $ 2 , 
ei(V$ ■ n) 1 - e 2 (V<I> • n) 2 = cr. 

En el caso de un conductor, dado que en su interior el campo electrico es 
nulo, se tiene 

^"conductor — C 011 ^? 
y asf, la densidad superficial de carga en el mismo es 

a = -e(V$ • n) afuera . 



Magnetostatica 

Las condiciones para hablar de magnetostatica son: 

• No hay dependencia temporal en los campos. 

• Hasta hoy no se han detectado los monopolos magneticos. 

Bajo est as consideraciones las ecuaciones de Maxwell ([l] - |J) se simplifican a 

VxH = j, (11) 
V-S = 0, (12) 

y las condiciones de frontera (|5| - H) 

{D 1 - D a ) • na, 2 = 0, 
n h2 x (H 2 - Hi) = z. 

Al igual que en el caso electrostatico, a partir del calculo vectorial y (|L?]) 
se introduce el potencial vectorial magnetico como 

5 = Vxl 

el cual, para materiales homogeneos e isotropicos (B = fiH), se obtiene de 
( |TTD como 

\/ 2 A = -fij (13) 

(NOTA: cabe aclarar que, dada su definicion, da lo mismo tomar A que 
A+V<p; por ello se elige el potencial vectorial tal que V(V- A) = 0, obteniendo 
asi (0) a partir de ([□]))• Si se conoce j, la solucion a (|13|) es 

Ait Jv \r — r | 

y para r > r sistema 

A ( r ) = 3 — 

donde m es el momento magnetico del sistema, dado como 

rh = - J rxj(r)dV. 
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Por ultimo, la energfa de un campo magnetico estatico es 
Wm*e = \ I B ■ HdV 



2 



mag 

f l{r) ■ . 
8n J |r — r*| 

Para un sistema de conductores 

Wmag = 2 X! L ik IiI k , 

i,k 

donde se define el coeficiente de induction magnetica entre las corrientes j$ 
y jfc como 

H 



Lik — 



4W, 



4J Ffc-rJ 
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Clase 2 



En el caso en el cual los campos varian lentamente en el tiempo, o sea son 
funciones f(at) donde se satisfacen las condiciones 

a < 1 uj < — I < A 

e 

con a c conductividad, uj, A caracteristicas de las oscilaciones electromagneticas, 

/ dimensiones lineales del sistema. 

Las ecuaciones de Maxwell toman la forma 

yxH = J V xE = —Qj: 

V-B = V'D = P- 

Observese que se ha despreciado el termino 

En el caso de campos variables arbitrarios para situaciones en las cuales no 
existen corrientes ni cargas presentes las ecuaciones de Maxwell toman la 
siguiente forma 

_ _ OB ^ dB 

vxH = ^ vxE = ~aF 

V-B = V'D = 0. 

Se tienen soluciones tipo ondas planas 

E = E e i(k ' r ^' ) 

H = H e i(k - r - hrt) . 



Las notaciones usadas son las usuales, uj es la frecuencia, I k 1= - es el 

1 'lie 

vector de onda, la direccion del cual, en medios isotropicos, coincide con la 
direccion de la energia. El vector que justamente nos da el flujo de energia 
es el llamado vector de Poynting (S), definido por 

S = ExH. (1) 

Para campos variables la conexion entre los campos y los potenciales es de 
la forma 
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B = v x A . (2) 

En general, los potenciales no son observables directamente (sino por sus 
efectos, E,B). Entre ellos existe una condicion muy importante (de consis- 
tencia de la teoria electromagnetica) que se llama condicion de "gauge" que 
puede ser diferente en funcion de la condicion considerada. 
Una de las condiciones de "gauge" mas frecuentes es la de Lorentz 

V'A + 6^=0. (3) 

Esta condicion "gauge" es muy usada porque permite una simple general- 
ization de las ecuaciones laplacianas del caso estatico 

□0 = -^ 
e 

□A = -fiJ (4) 

estas ecuaciones son llamadas D'Alembertianas y □ = y 2 — J^-. Los po- 
tenciales que son solution de estas ecuaciones son llamados potenciales re- 
tardados (y no precisamente porque sean muy tontos) los cuales tienen la 
forma 

A{r ' i) = il 3W \l-X ) dV '' (6) 

A grandes distancias del sistema de cargas (r >> A) y en el vatio, B,E y A 
se pueden escribir como sigue 

B = -Axn (7) 

c 

E = cB x n = (A x n) x n (8) 

A = / J(r',t- 1 r ~ r/ \- )dV' (9) 
Anr J v 

donde n= - es el versor en la direccion de la radiacion. Si ademas A » L con 

r 

/ la dimension del sistema radiante, se puede usar la llamada aproximacion 
multipolar, es decir, la radiacion se puede representar como una sumatoria de 
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los campos emitidos por los dipolos, cuadrupolos, etc., que forman el sistema. 
Para el caso dipolar se tiene 



^ P x n , . 

B = "°l^ < 10 > 

E=£-(px„x„) (11) 

donde p es el momento dipolar del sistema. La intensidad de la radiacion de 
un dipolo es 

' = «ife- (12) 

Magnetohidrodinamica 

La magnetohidrodinamica estudia el comportamiento de los liquidos o los 
gases conductores (plasmas) en campos electromagneticos. Se usan los con- 
ceptos hidrodinamicos: densidad, velocidad, presion, viscosidad. Las ecua- 
ciones basicas son: 

d P m , „ ( \ n 
-^T + V • M = 

dv , 2 

Pm-r^ + Pm(v ' V) V = -V^+J xB + ^V V + Pmg 

yxE = -— yxH=j j = u e (E + vxB) 
mas la ecuacion de estado del fluido. 



Relatividad Especial 

La teoria de la relatividad especial surgio en la electrodinamica y se basa en 
dos postulados fundamentales 

• La velocidad de la hrz en el vacio c = 2.99793 x 10 8 m/s es una constante 
en todos los sistemas de referenda inerciales. 

• Las leyes de la Fisica tienen la misma forma en todos los sistemas 
inerciales (covariancia de las leyes naturales). 
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Las transformaciones de Lorentz en una dimension se escriben asi 



x\ + if3x± 



X2 — X2 



X A 



X4 



X 3 — £3 

if3x\ 



donde x 4 = ict, x' A = ict', (5 — -. Las velocidades vl de un cuerpo en K' 
con respecto a las velocidades u del mismo cuerpo en K estan relacionadas 
mediante las siguientes expresiones 



u., 



^ _ VUx ' 



Uy = 



Y vu x ' 



u z = 



La segunda ley para particulas relativistas se escribe 



dp 

~dt 



d 



mv 



dt \VT=P 



Cantidades del tipo (p,i/cE) son llamadas cuadrivectores, el anterior se 
llama cuadrivector de impulso-energia. Otros ejemplos de cuadrivectores son 
(k, i/cu); (j,icp); (A,i/cip). Existen tambien objetos llamados cuadriten- 
sores por extension de lo anterior, algunos ejemplos de ellos son 



a/3 



( 





cB z 


-CBy 


-iE x \ 




-cB z 





cB x 


-iE y 




cB y 


-cB x 





-iE z 


V 


%E X 


lEy 


iE z 


/ 



T a /3 — ^o(F af iFf3n — -^apF^F 1 * 71 ) 
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Clase 3 



Fuerza de Lorentz como fuerza lagrangiana 

Las ecuaciones del movimiento de Euler-Lagrange son 

Q ~ -— + — (—) (1) 

k dq k dt \dq k ) 

donde las Q k son las fuerzas externas o fuerzas generalizadas y L = T — U . 
Por otra parte, las ecuaciones de Maxwell en unidades de Gauss son 

(Ml) V x ~E + = (MS) V-L)=Anp 

IdD An—* 

(M2) VxiJ — = — j (MA) V • B = . 

C CJ v c 

Ahora con F = q E = —qV(f solo en electrostatica en general la fuerza 
es la ley de Lorentz o sea 

Fi^p + ^xB) . (2) 

Ahora de la (MA) encontramos que B = V x A y sustituyendo en (Ml) 
encontramos 

V x ~E +-— (V x ~A ) =0 (3) 



cdt 
por tanto 

de aqui que podemos definir una funcion escalar tal que 

-V$ = ^ + --A (5) 

c at 

entonces 

t l = q U*-\p + Wxhx^)) (6) 
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donde el doble producto vectorial lo podemos expresar de la siguiente forma 
V x ( V x 4 ] = V ( T ■ A ) - + — (7) 



dt dt 



por tanto 

~Fl = q 

= g(-v 



V 



COT C 



+ 



$ V • A 

c 



1 d 

c rft 



dt eft 
• A 



lo que hace que F l se pueda escribir como fuerza lagrangiana 

d dU 



F l — -VU + 



dt Q~ 



(8) 



con U — q® — 1 ~v ■ A . 



Electrodinamica no lineal 

Para la electrodinamica no lineal la constante dielectrica se expresa como 

e v = t — — rv (9) 

y la permeabilidad se escribe como 

^ = ^o{l + ^ 2 (c 2 B 2 -E^ (10) 

donde b en ambos casos es un parametro que fija una intensidad maxima de 
los campos. 

Al menos para campos que varfan lentamente, en funcion de los tensores 
de permeabilidades electricas y magnetica del vacio tenemos 

Di = ^2e ik E k y Bi = ^2fi ik H k (11) 

k k 
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donde 



&ik 



457rm 4 c 



-2 (E 2 - c 2 B 2 ) 5 ik + 7c 2 B t B k 



Ahum^c 7 



2 B 2 



E 2 



5 ik + 7EiE k /c 2 



+ ... 
+ .. 



(12) 
(13) 



para el lmiite clasico hacemos h ^ j estos efectos no lineales desaparecen 
al comparar con la expresion clasica en (9) y (10) encontramos 



0.51- 



por tanto 



V457T / e 2 e 
2 y A-KEohcA-KSorl " ' AKE r1 



e 2 

r a = -2- « 2.8 x 10~ 15 metros 
mc 2 



0.51 



(14) 



(15) 



este es el radio clasico del electron. 
Ahora si tenemos varias cargas 



E -~n da = — ^ ft 



(16) 



y si tenemos distribuciones de carga 



1 



E ■ n da = — / p (x ) dV 

£ JV 



(17) 



donde V es el volumen enserrado por la superficie, ahora el teorema de la 
divergencia nos dice que 



v ■ n da 



V • v dV 



v 



entonces aplicando este teorema en la ley de Gauss encontramos 

V-E=P- 



(18) 



(19) 



y esta es la forma diferencial de la ley de Gauss. 
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Clase 5 



Energfa potencial electrostatica y densidad de energfa; 
capacitancia 

Imaginemos el caso en que una carga % es traida desde al infinite hasta el 
punto Xi, localizado en una region del espacio donde se conoce el potencial 
electrostatico El trabajo realizado sobre esta carga es 

Wi = qMxi). 

Ahora bien, si este potencial es provocado por la presencia de otras n — 1 
car gas, se tiene 

1 n_1 a 



47re ~[ \Xi-Xj\ 

j por tanto 

^i^E^ (!) 
47re o ^ Fi - Xj\ 

Por un proceso mental similar, se puede ver que el trabajo total necesario 
para obtener el arreglo de n cargas, trayendo cada una desde infinite a una 
region del espacio originalmente vacia, es 

donde i,j toman todos los valores entre lyn, excepto i — j (autoenergias). 
En el caso de una distribucion continua de cargas es claro que 

expresion que puede reescribirse de varias formas: 
• En terminos del potencial 

Wtotal = \ j P {xmx)<fx. (4) 
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Utilizando la ecuacion de Poisson: 



Wtotal = "I / ®Wx. (5) 

Integrando por partes la ultima expresion se obtiene 

Wtotal = |/|V$| 2 rf 3 x 

= || |£| 2 d 3 x. (6) 

Por la forma de la ultima integral, se define la densidad volumetrica de energia 
como 

Notemos que esta densidad de energia es no negativa, y por tanto el tra- 
bajo total tampoco sera negative Sin embargo, de (|J) se ve que el trabajo 
para hacer un arreglo con dos cargas de signo contrario es negativo; esta con- 
tradiction surge porque en las expresiones (|3]- ^) se incluyen las autoenergias 
en el trabajo total, mientras que en el caso discreto (g) se las excluye. 

Por ultimo, como siempre, se puede calcular la fuerza a partir de los 
cambios que sufre la energia ante desplazamientos virtuales pequefios. 

Consideremos un sistema de n conductores, el %— esimo de ellos con carga 
Qi y potencial Vi. Dada la relacion lineal que existe entre el potencial y la 
carga, podemos escribir 

n 

i=i 

donde Pij depende solo del arreglo geometrico de los conductores. Invirtiendo 
las ecuaciones anteriores se obtiene 

n 

Qj = ^2 CjiVi- 
i=l 

Los coeficientes Ca son las capacitancias, y (i ^ j) los coeficientes de 
induccion. 



De esta forma 



^total = \flQiVi 
z i=i 
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Aproximacion variacional a la solucion de las ecua- 
ciones de Poisson y Laplace 

El uso de metodos variacionales es muy popular en Ffsica. La electrodinamica 
no es la exception. En efecto, la idea de considerar funcionales cuyos ex- 
tremales satisfagan ecuaciones de movimiento tipo Poisson o Laplace es muy 
sugestiva (sobre todo por la elegancia del metodo variacional). 
Consideremos la funcional 

Ity] = j VV> • Vtpd 3 x - f gipd 3 x, (7) 

sujeta a la condition tipo Dirichlet 6if)(S) = (S es la superficie cerrada que 
contiene a V). Es facil ver que 51 = + Sip] — = conduce a la 
ecuacion de movimiento 

V 2 ^ = -g. 

Se ve que este problema no es otro que resolver la ecuacion de Poisson con 
condiciones de frontera tipo Dirichlet. 

Similarmente, para condiciones de frontera tipo Neumann, se plantea el 
funcional 

IU] = - f • V4>d 3 x - f gipd 3 x - I fipd 3 x, (8) 
2 Jv Jv Js 

con 

Es facil probar que SI[ip] = conduce a las ecuaciones 

V 2 t£. = —g, 



dn 



s 



Resulta logico preguntar si este metodo variacional de obtener la ecuacion 
de Poisson sirve para algo, o es solo un juego matematico. Para contestar, 
notemos que una vez conocida la forma de los funcionales (0, ||) aiin es 
necesario encontrar ip (o sea resolver la ecuacion de Poisson); por tanto el 
problema es el mismo. Sin embargo, se pueden proponer soluciones ip = 
^^(a;, a, j3, ...) que satisfagan las condiciones de frontera dadas, para despues 
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variacionalmente encontrar las constantes indeterminadas (notar que con esta 
election I = I[A,a,(3, ...]). En este sentido el metodo variacional sirve para 
encontrar soluciones aproximadas. 
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Clase 6 



Metodo de las imageries 

Este metodo se refiere a problemas de cargas puntuales en la presencia de 
superficies a potencial cero o constante. Las condiciones de frontera se simu- 
lan con cargas puntuales de valores y posiciones bien determinadas conocidas 
como "cargas imagenes". 



Carga puntual con esfera a = 

El potencial asociado a la carga real y la carga imagen es 



47re | x - y | | x - y' | 

La condition de frontera es que el potencial se anule en | x |= a. Introducimos 
dos vectores unitarios n,n', uno en la direccion de x y el otro en la direccion 
de y, de manera que el potencial se puede expresar 

a( \ 1 r q i i 



47re | xn — yri \ \ xn — yn 1 \ 
Factorizando x del primer termino, y' del segundo y valuando en x = a 

fa = a ) = —I— TTT-\ + 77TZT/ ^ 



47re a \ n — jn' \ y' | n' — ^-n 



Se observa que para que el potencial se anule en la frontera de la esfera se 
debe satisfacer 

q _ _4_ y_ = ^_ 

a y' ' a y' 

resolviendo estas ecuaciones se encuentra 

, y' a , a 2 

q = q = q, y = — 

ay y 
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q' es la carga total de induction sobre la superficie de la esfera, podemos 
observar ademas lo siguiente 

y -> a q' -> -q 
y — > oo =>- </ — > 
La densidad superficial de carga esta dada por 



a = -e 



g a 



y 1 



dx 



47ra 2 y(l + ^-2^cos7) 3 / 2 

Es posible calcular tambien la fuerza de atraccion hacia la esfera, la magnitud 
de la cual esta dada por 



2 „3 



1 q 2 a 
47re a 2 y 3 



1 - 



Carga q en presencia de una esfera conductora cargada 
a Q, aislada 

El potential para esta configuration se puede expresar asf 



0(x) 



4vre 



aq 



x-y 



+ 



Q 



y I x 



La fuerza de atraccion en este caso es 



F(y) = 



i gy 

47re y 3 



Q 



qa 3 (2y 2 - a 2 ) 
y(y 2 - a 2 ) 2 



Carga q cerca de una esfera conductora a potencial 
constante 



Para la situation presente el potencial adopta la forma 



0(x) 



47re n 



aq 



x-y 



y I x - ^y 



+ 



Va 



La fuerza de atraccion esta dada por 



F(y) 



<?y 



Va- 



qay° 



47re (y 2 - a 2 y 
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Esfera conductora en un campo electrico uniforme 

Un campo electrico uniforme es producido por ejemplo por dos cargas pun- 
tuales ±Q localizadas en z = ±R para R — > oo. Si ahora una esfera conduc- 
tora es colocada en el origen, el potencial sera el debido a las cargas ±Q en 
q=R y sus imagenes =F^ en z = =F^ 



= 



Q 



Q 



+ R 2 + 2rR COS0) 1 / 2 (r 2 + R 2 - 2rR cos 6) x l 2 



+ 



47re r 



aQ 



+ 



aQ 



R(r 2 + ^ + 2|r C os^) i?( r 2 + ^_ 2|r cos ^ 

Como R» r podemos desarrollar los denominadores 

1 



1/2 







47re 



2Q . 2Qa 3 
-— r cos 6^ + ^ cos ft 
it 2 it 2 r 2 



+ 



Para it! — > oo, 4n ^ R2 es el campo aplicado de manera que el potencial en este 
Kmite toma la forma 

/ a 3 \ a 3 
(f)R^oo = -E ( r - — J cos 9 = -E z + ^E z , 

donde el ultimo termino es el del "dipolo imagen". La densidad superficial 
de carga esta dada por 



a = -e — \ r =a= 3e E cos9 , 
or 



la cual se anula al integrarla sobre la superficie 

J ada = . 
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Funcion de Green para la esfera conductora 

Para problemas de Dirichlet con conductores G(x, x')/47reo puede ser inter- 
pretada como el potencial debido a la distribucion superficial de carga in- 
ducida sobre la superficie por la presencia de una carga puntual (fuente) en 
el punto x'. Por definition la funcion de Green G(x, x') satisface la ecuacion 

V /2 G , (x,x') = -4tt5(x-x') . 

Para el caso de la esfera la funcion de Green esta dada por 

1 a 



GW(x,x') = 



En coordenadas esfericas lo anterior es 



rtf' -\r- ~ ~V~' 

I x l2 A 



G ' 1 1 

(x 2 + x' 2 -2xx'cos-f) 1/2 ^ + a 2 -2xa;'cos7) 1/2 



dG , ~ 2 - ~ 2 



x — a 



dn' a (x 2 + a 2 — 2ax cos 7) 3//2 

Recordando la solucion de la ecuacion de Poisson con condiciones de Dirichlet 
para el potencial 

1 r If r)G 

usando esto, podemos escribir la solucion general para el potencial de la 
esfera conductora para la cual conocemos el potencial en la frontera 

4vr j ( x 2 + a 2 _ 2ax cos 7) 7 

donde cos 7 = cos 9 cos 9' + sin sin cos(</? — <//). Para el interior de la 
esfera x 2 — a 2 — > a 2 — x 2 , y en el caso en el que se tienen distribuciones 
volumetricas de carga se tiene que tomar en cuenta la contribution de la 
integral de volumen. 
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Clase 9 



Analisis de elemento finito para resolver la ecuacion 

de Poisson 

A continuation presentamos una breve introduction al analisis de ele- 
mento finito para resolver la ecuacion de Poisson. Por simplicidad en la 
presentation solo consideramos problemas bidimensionales. 

Primeramente esbozamos el metodo de Galerkin para replantear la ecua- 
cion de Poisson, y dividir la region de estudio en una red cuyo numero de 
celdas es finito. Por ultimo presentamos dos tipos particulares de redes: 
cuadriculada regular y triangular. 

El metodo de Galerkin 

Sea una region bidimensional R limitada por una curva cerrada C; consider- 
emos en R la ecuacion de Poisson 

V 2 V = -9 (1) 

con condiciones de frontera tipo Dirichlet; multiplicamos (HD por una funcion 
de prueba <f>(x, y) que sea continua a trozos en R y tal que 4>{C) = 0; despues 
integramos sobre R, obteniendo 

/ [0V 2 ^ + g4>]dxdy = 0. 

JR 

A continucion, utilizando la primera identidad de Green (bidimensional), la 
integral anterior se reescribe como 

/ [V0 • - g<p}dxdy = 0. (2) 

•J R 

El siguiente paso es dividir la region R por medio de una red con iV celdas, 
y definir un conjunto de funciones {<fii(x, y),i = 1,2,..., iV} tal que cada una 
de ellas es no nula solo en una celda particular de la red. A continuation se 
expresa ip como 

N 

^(x,y) w S*i0i(ar,j/); 
i=i 
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sustituyendo lo anterior en (0) y escogiendo = <pj se obtiene 

/ V0i(ar,j/) • V(/)j(x,y) = go 4>i(x,y)dxdy , 
i=i J Ji? 

donde se ha supuesto que las celdas son suficientemente pequenas como para 
que g(x,y) ~ g dentro de ellas (el valor de g varia de celda a celda). Con 
esto, @) se reduce a la ecuacion matricial 

K\l> = G (3) 

aqui K es una matriz N x N con elementos 

kij = / V0i • V jdxdy 

J Ft 

\l/ es la matriz columna formada con los coeficientes ^f, G es una matriz 
columna con elementos 



Gi = gi / <f)i{x,y)dxdy. 
Jr 



El poder del metodo de Galerkin radica en que, por la forma como se escojen 
las (pi, la matriz K es dispersa, i.e., solo pocos de sus elementos son diferentes 
de cero, y por ello es relativamente facil conocer \1/ a partir de (||), lo cual 
nos da la solucion a la ecuacion tipo Poisson ([[]). 



Casos particulares 

Red cuadriculada regular 

Se escoge una red de cuadros, cada uno de lado h; sean (xj, yj) las coordenadas 
de los vertices. Se toman las funciones <pij(x,y) tales que 0^ ^ solo en 
una vecindad de area h 2 alrededor de (xi,yj), y las son linealmente 
independientes entre si. Con esto, de acuerdo al metodo de Galerkin 

(No) 

ip « ®ki(t>ki(x,y) 

k,l=l 
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donde se supone que el total de celdas es N ; los coeficientes tyf-i se obtienen 
a partir de @ con 



K = ^ V0y ■ Vcpudxdy 
(G) = (g{x h yj) I fajdxdy 



(*) = (*<)■ 

La inconveniencia del uso de redes como est a es que se presentan casos 
donde el potencial varia de formas diferentes en diferentes regiones, y por ello 
seria mas conveniente utilizar celdas irregulares. A continuation se presenta 
una de ellas. 



Red triangular 

Las redes triangulares son las mas utilizadas en el analisis de elemento finito, 
por las razones expuestas al final de la section anterior. Para este tipo de 
redes se asume que el elemento triangular (e) es lo suficientemente pequeiio 
como para que ip cambie poco en su interior y de hecho pueda ser aproximado 
de forma lineal en cada direction: 

if)(x, y) ~ if> e ( x , y) — A + Bx + Cy. 

Sean (xj,^) (i = 1,2,3) las coordenadas de cada vertice del triangulo. En- 
tonces las constantes (A, B, C) quedan determinadas por los valores de if) en 
cada uno de ellos. 

Con el fin de sistematizar el procedimiento, es conveniente definir las 
funciones de forma Nj(x,y) (una por cada vertice), tales que Nj e \xj, yj) = 1, 

si x 7^ Xj, y 7^ yj. Por la linealidad de if) 

dentro de e, tomamos Nj(x, y) = ctj + bjx + Cjy. De aqui, para j — 1 

ai + biXi + c x yi = 1 

ai + bix 2 + c\y 2 = (4) 
«i + hx 3 + cij/3 = 

de donde 

a i = 7^( X 2V3 ~ x 3 y 2 ) 
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donde S e es el del triangulo e. 

Ahora, siguiendo el metodo de Galerkin, tomamos fa = Nj; e \ De esta 
forma, expresamos ip como 

^^EW^j/). ( 5 ) 

donde la suma se realiza para todos los triangulos (/) y todos los vertices de 
cada triangulo (j); \E^- es el valor de ip en el vertice j del triangulo /. Estos 
coeficientes se encuentran, para cada triangulo, a partir de una ecuacion 
similar a @: 

Eff=k^ 
i=i 6 

con fcjj = S e (bibj + CiCj) (coeficentes de acoplamiento) ; g e = g(x e ,y e ), y 
(x e , y e ) son las coordenadas del centro de gravedad del triangulo. A contin- 
uation solo falta incluir todos los triangulos de la red. Para ello, considerando 
que los vertices interiores a C son N, y el total de vertices (interiores a C 
y sobre ella) es N Q , los indices corren de 1 a iV para los vertices internos, y 
de N + 1 a N para los que estan sobre la frontera. Con esto, se obtiene la 
ecuacion equivalente a fl3|) para toda la red es con 

K = (kij), ku = E ■> kij = E kij, i 7^ j, 

T E 

I ^0 

^ = ^E^ e - E *#*; e) ; 

d T J'=JV+1 

T indica que la suma es sobre los triangulos con vertice comun %\ E que la 
suma es sobre triangulos con lados entre los vertices i, j. 

Como ya se dijo, K es una matriz dispersa, y por tanto la solution a ([I]) 
se puede encontrar como 

f>3 
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